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q \ It turns out that in any case the purported relation to differentially algebraic 

functions, and hence to Shannon's model of analog computation, fails. 

1. Introduction 

m 

There are several different kinds of theoretical models that talk about "computability" and 



Moore introduced a class of real-valued "recursive" functions by analogy with 
Kleene's formulation of the standard recursive functions. While his concise def- 
inition inspired a new line of research on analog computation, it contains some 
technical inaccuracies. Focusing on his "primitive recursive" functions, we pin 
down what is problematic and discuss possible attempts to remove the ambiguity 
regarding the behavior of the differential recursion operator on partial functions. 



"complexity" of real functions. Computable Analysis 19( and some other equivalent models 
use approximation in one way or another to bring real numbers into the framework of the 
standard Computability Theory that deals with discrete data in discrete time. Another well- 
known model is the Blum-Shub-Smale model [lj] in which continuous quantities are treated 
as an entity in themselves but the machine still works with discrete clock ticks. 

A third approach is analog computation in which not only are the data real-valued, but 



also the transition takes place in continuous time [15|. One of the oldest and the best- 
studied models of such computation is Shannon's General Purpose Analog Computer 18 
that models the Differential Analyzer [sj, a computing device built anclput to use during 
the thirties through the fifties. The GPAC, after some refinements jsl, H2, 16], was shown 



capable of generating (in a sense) all and only the differentially algebraic functions. We will 
explore this class in Section [5] and show that it can be characterized in many different ways. 
Little is known about how such analog models relate to the standard (digital) computabil- 



ity. Moore [13|] addressed this question for his new function classes that also try to express 
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the power of GPAC-like computation. In imitation of Kleene's characterization of the usual 
recursive functions, these classes are defined as the closures under certain operators that 
are supposedly real-number versions of primitive recursion and minimization. He makes the 
following claims, among others, that relate his classes of real primitive recursive and real re- 



cursive functions to analog and digital computation, respectively |13|, Propositions 9 and 13]. 



Claim 1. Real primitive recursive functions are differentially algebraic^ 

Claim 2. Each (partial) recursive function on the nonnegative integers (in the standard 
sense) is a restriction of some real recursive function. 

Despite its impact on the subsequent study on the classes and their variants 
his work lacked formality in some ways, as already pointed out 0,0]. In fact, the definition of 
the classes suffers from ambiguity. In Section [3] of this paper, we reformulate Moore's theory 
up to the real primitive recursive functions in a mathematically sound way. With the aid of 
the preparation in Section El we show that, even though our formulation of the class seems 
the most restrictive possible, Claim CD fails. In section HI we discuss some issues about classes 
other than the real primitive recursive functions, including Claim [2j 

Throughout the paper, we write N, Z, R for the sets of nonnegative integers (including 
0), integers and real numbers, respectively. 

Partial functions In this paper, a function f :C R m — > R n may be partial, as opposed to 
total; that is, the set dom/ of x G R m for which the value f x & R" is defined is allowed to 
be a proper subset of R m . By the restriction of / to a set J C R m we mean the function g 
with doing = J fl dom/ such that gx = fx for every x G domg. When dom/ is open, / 
is said to be (real) analytic if for every a = (a , . . . ,a m -i) £ dom/ there are an open set 
J C dom / containing a and a family (c p ) pg N m of ^-tuples of real numbers such that the sum 

c p ■ (x - a ) pa ■ ■ ■ (x m _i - a m ^) Pm - 1 (1) 

P=(p ,---,Pm-l)eN m 

converges to fx for each x = (xq, ■ ■ ■ ,x m -i) G J (regardless of the ordering of summation). 
See Krantz and Parks (lO . Chapters 1 and 2] for well-known properties of analytic functions. 
When / is analytic, we write d^ 00 '™' "- 1 )/ (and not d ao+ - +am - 1 f /dt%° ■ ■ ■ dt^) for the 
mixed partial derivative of / of order aj along the i th place (which is known to exist). 



Moore [13] does not explicitly deal with partial functions. We believe that this is responsible 
for ambiguous and erroneous statements made in his seminal work as well as in some of the 
subsequent works by other authors. Although there are some situations in mathematical 
analysis where we can pretend that there are no partial functions (namely, when we are only 
discussing properties defined locally, such as continuity or analyticity), this is not the case 
with the notions we want to discuss here. If, say, the above Claim [5] is to make any nontrivial 
sense, it is clearly inappropriate to talk about "real recursiveness at x," as there is a real 
function which is simple locally but the restriction of which to N is highly complicated in the 
recursion-theoretic sense. We therefore emphasize that partial functions must be dealt with 
seriously, and devote this paper to accordingly reformulating the theory wherever possible. 



^^Moore writes Mq for the class of real primitive recursion functions. Claim[T]was later replaced by a similar 
claim [1, Proposition 2] for a more "restricted" class Q than M , but its definition is again unclear. 
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2. Differentially algebraic functions 



We show some facts about single- and multi-place differentially algebraic functions. 

Theorem 3. Let m, n and i < m be nonnegative integers and f :C R m — > R n be an analytic 
function with open domain. Let ®, (Jn]), ([m|) and ffrvT) be the following statements: 

(i) for any open connected set J C dom/, there is a Z™ - coefficient nonzero polynomial P 
such that 

P (fx, D e </x, D 2ei fx, B {aT[tyP - 1] ^fx) = (2) 
for all x £ J, where G N m zs ine vector whose i th component is 1 and others are 0; 

(ii) for each Xq G dom /, inere are a Z n - coefficient nonzero polynomial P and an open set J 
containing xq such that we have (EJ) for all x G J; 

(hi) for each xq G dom/, tnere are a Z" - coefficient nonzero polynomial P and an open 
interval J containing the i th component of xq such that we have (T5]) for all x whose i th 
component is in J and whose other components equal those of x$; 

(iv) for each x G dom/, there is a Z™ -coefficient nonzero polynomial P satisfying (|2"1) . 

Let (Eh), (Elk) and (lmk) fre tae statements obtained by replacing 7i by R in (pj, (El) ana 1 (lull) . 
respectively. Then (0), (jn]), dm]), (jrvj) . (Jib); (JUtt) (jSIk) are equivalent. 

Proof. The implications (0) ^ (Ell) ^ (En]) =>• (Evj) and (EJ) ^> (E^,) ^> (JTTtt.) =>■ (EHk) are obvious. 
It has been known that (lmk) =>• (pHl) . see Theorem [TBI To see (ITv]) ^> ([!]), consider, for each Z n - 
coefficient polynomial P, the set Jp of all x G J satisfying ([2]). Since by (JIvJ) these countably 
many closed sets Jp cover the open set J, one of them must have nonempty interior by Baire 
Category Theorem [TH This Jp must then equal J by the Identity Theorem [15j □ 

Definition 4. Let m and n be nonnegative integers. An analytic function / :C R m — > R n is 
differentially algebraic^ if for each i < m it satisfies one (or all) of the clauses in Theorem [3j 

Note that / need not be the unique solution of (EJ). For example, every function (with 
open domain) that is constant on each connected component of its domain is differentially 
algebraic because of the single set of equations D £i fx = 0. 

The clauses (pi])- (pvj) show that being differentially algebraic is a "local" property. 

When dom / is connected, @ reduces to the following statement: 

d3) there is a Z n -coefficient nonzero polynomial P such that we have ([2]) for all x G dom/. 

Hence, the clause (En]) shows that, as long as dom / is connected, our definition is equivalent 
to that of many authors, including Moore 13J, who first state the definition for m = 1 by ([!]) 
and then extend it to m > 1 by saying that a function is differentially algebraic when it is so 
as a unary function of each argument when all other arguments are held fixedJl We proved 
Theorem [3] in order to use (EJ) to present a counterexample to Claim [TJ later. 



2 Also termed algebraic transcendental or hypotranscendental. Functions without this property is said to be 

trans cendentally transcendental or hypertranscendental. 
3 Their definition for the case m — 1 is slightly different from (jU) in that it replaces ^ by P(x, fx, D ei fx, . . . , 

D^ arlty p ~ 2 ^' ei fx) — 0. But the proof of (jaj => JbJ of Lemma [5] shows that this difference is superficial. 
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Let us characterize differentially algebraic functions in yet another way for the case n—1. 
For a field E, its subfield F and a set B C E, we write F(B), agreeing tacitly on E, for the 
smallest subfield of E that includes F and B. We write F for the algebraic closure of F, that 
is, the set of those elements of E that annuls some F-coefncient unary nonzero polynomial. 

Let J C R m be an open set and consider the ring C W [J] of analytic functions g :C R m — > R 
with doing = J. Note that R is embedded into this ring by regarding each x G R as the 
constant function taking the value x. To assert (J2J) for all x G J is to say that 

P (f, D e 7, D 2 ' e 7, . . . , v(^yP-V^f) = o (3) 

in C w [J] . If J is connected, C u [J] has a quotient field by the Identity Theorem [151 so the 
notation R(D /) in the following lemma makes sense. We write D / = { D a f | a G N m }. 

Lemma 5. Let J C R m be open and connected. For f G C W [J], the following are equivalent: 

(a) / is differentially algebraic; 

(b) D / C R(B) for some finite set B CD f; 



(c) D / C R(B) for some finite set B C R(D /). 

Proof. The implication (jbj) =>■ (jcj) is trivial. The Transcendence Degree Theorem [I7J shows 
(Ej) =>■ ((Ej). For ([aj) =>• (jbj), assume that for each i we have an R-coefficient polynomial Pi 
with 

P f (/,D«/,D 2 ^,...,D^/) = 0, (4) 

where iVj = arity Pj — 1 . By choosing iVj to be smallest and then the degree of Pi in the last 
place to be smallest, we may assume that 

S = ( D (0,...,0,l)p. ) (/> D e l/; D 2- ei/) ^ D iV, ei/) (5) 

is nonzero. Consider the order < on N m defined by setting a < b when a + c = b for some 
c G N m . We will show by induction on a G N m that D a / G R({ B b f \ b< (JV , • • • , A^ m _!) }). 
The case a < (A/q, . . . , N m _i) being trivial, assume that a > (Ni + 1) ■ ej for some i. Apply 
Y^a-Ni-ei ^ Q Ytoth. sides of (jlj) and calculate using chain rules to obtain 

& + E -D a f = 0, (6) 

where \P can be written as a sum of products of several derivatives of / of order < a and 
a, which hence enjoy the induction hypothesis. □ 

Apart from purely theoretical interest, the significance of differentially algebraic functions 
lies in their relation to the General Purpose Analog Computer, an analog computation model 



introduced by Shannon [18( and later refined by Pour-El [16(. More precisely, if a function 
/ :C R — ► R with nonempty domain is differentially algebraic, then the restriction of / to 
some nonempty subset of dom / is GPAC generable [la Theorem 4] ; conversely, if a function 
/ :C R — > R with nonempty domain is GPAC generable, then the restriction of / to some 
nonempty subset of dom / is differentially algebraic 0, Theorem 2] . Graga later considered 
the Polynomial GPAC, a simpler refinement than Pour-El's, and proved analogous results jsj. 
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3. Real primitive recursive functions 



The class of real primitive recursive functions is denned [13j as the smallest class containing 
some basic functions and closed under the operators specified below. Unfortunately, the 
original definition contains ambiguity, resulting in some inconsistent claims about the class. 
To remedy this, we shall revisit the definitions carefully in Sections 13.11 and 13.21 Section 13.31 
discusses an alternative approach by Campagnolo. Section 13.41 disproves Claim [TJ 



3.1. Two basic operators 

The real primitive recursive functions are defined through three operators: juxtaposition, 
composition and differential recursion. The first two are very simple. 

Definition 6. Given go, ... , g n -\ :C R m — > R, define their juxtaposition JX (go, . . . , g n -i) 'Q 
R m — > R n by setting JX (g , . . . , g n -i)x = (gox, . . . , g n -ix) whenever x G dom g R - • • H 
dom(7 n _!. Given / :C R m — > R n and g :C R' — ► R m , define their composition CM (/, g) :C 
R' — ► R n by setting CM (f,g)x = f(gx) whenever x G doing and gx G dom/. 

We write fog for CM (f,g). As remarked in Section [U it is important to define precisely 
what the operators do on partial functions. Note how Definition [6] specifies the domain of 
the functions constructed. If gx is not defined, neither is (fog)x, even if, say, / is a constant 
function defined everywhere. We thus work in the following (informal) general principle. 

Principle 7. For the value of an expression to be defined, the value of each of its subexpres- 
sions has to be defined. 

We remark that this was not explicitly intended by Moore. In fact, he presents an example 
to the contrary when he claims 13|, Section 6] that the total function inv :C R — > R given 
by 

— JO ifx = 

mvx — < (7) 
[l/x ifx^O W 

can be obtained by composing the binary multiplication with JX (zero?, g), where zero? is 
from (IHUjl and g is the restriction of inv to R\{0}. Some authors point this out [1, p. 22] and 
criticize it [7], p. 47]. Without discussing which definition is more "natural," we adopt our 
restrictive Definition [6j simply because it is not clear how to formulate a general definition 
that would admit this construction of inv. 

The operators JX and CM preserve analyticity Proposition 2.2.8]. 

Theorem 8. The property of being differentially algebraic is preserved by JX and CM. 

Proof. This is trivial for JX. For CM, it suffices to show that if / :C R m — > R and go, ■ ■ ■ , 
flVn-i :C R' — > R are differentially algebraic, so is / o g, where g = JX (g , . . . ,g m -i ). We 
may assume that doing and J = dom (/ o g) are connected. We use the characterization (jb]) 
of Lemma EJ Calculate each element of D (/ o g) by the chain rule to see that it belongs to 

m—l 

K({dog\deBf}U \J{q\j\qeT> gi }), (8) 

where q \j means the restriction of q to J. By the assumption, there are finite subsets ACD/ 
and Bi C D gi with D / C R(A) and D C R(_Bj) for each % — 0, . . . , m — 1. This implies 
that (jHJ) stays unchanged by replacing D f by A and D ^ by □ 
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Figure 1: When the equation is satisfied by both h and hi (as well as their restriction to 
each interval containing the origin), how do we say that the shaded interval is where 
there is a "unique solution"? 

3.2. The differential recursion operator 

To formulate the third operator, we need a notion of unique solution of an integral equation 
of the form ([9]) below, where h is the unknown. For example, it sounds natural to say that 
the tangent function restricted to (— tt/2, 7r/2) uniquely solves ht = + (hr) 2 ) dr. But as 
we are talking about partial functions, the word "unique" should be used carefully, because 
the restriction to any subinterval J C (—tt/2, tt/2) containing also satisfies the equation on 
J. Thus, out of the set H of all solutions, we need to pick one function that deserves to be 
called the unique solution defined on the largest possible interval (Figured]). Though Moore 
did not discuss this, it is not hard to formulate this intuition: for a set H of functions of a 
type, we say that a function h G H is unique in H if the restriction of any function in H to 
dom h is a restriction of h. 

Definition 9. Let / :C R m -> R n and g :C R™+i+™ _> R™. For each v G R m , let H v be the 
set of all functions h :C R — > R n such that 

(a) dom/i is either the empty set or a possibly unbounded interval containing 0, 

(b) v G dom / if dom h is nonempty, 

(c) (v, r, hr) G doing for each r G dom/i, and 

(d) every t G dom h satisfies 

ht = fv+ / g(v,r,hr)dr. (9) 
Jo 

Let K v be the set of functions unique in H v . By Lemma [18] in the appendix, K v has an 
element h v of which all functions in K v is a restriction. Define DR (f,g) :C R m+1 — > R n by 
dom (dr (/, g)) = { (v, t) G R m+1 | t G dom h v } and DR (/, g) (v, t) = h v t. 

Definition 10. The class of real primitive recursive functions is the smallest class containing 
the miliary functions O ^ 1 , l ^ 1 , — l ^ 1 and closed under JX, CM and DR. 
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Lemma 11. The following functions are real primitive recursive: for each neN, the n-ary 
constants O™^ 1 , \ n ^ 1 , — l n ~* 1 ; for n G N and i — 0, . . . , n — 1, the n-ary projection id™ - * 1 to 
the i th component; binary add and mul; the functions inv + (mapping x > to 1/x), sqrt + 
(mapping x > to y/x) and In (natural logarithm) defined on (0, oo); the total functions sin, 
cos and exp; the circle ratio it as a nullary function. 

Proof. The constant O™ - * 1 is built by n->1 = O ^ 1 o jx(); similarly for l™^ 1 and — l n_>1 . 
Then inductively define id^ 1 ^ 1 = DR (0^\ l^ 2 "* 1 ) and id™ +1 ^ = DR (idp 1 , n+2 ^ 1 ). Using 
these, let add = DR^idjif* 1 , l 3 ^ 1 ) and mul = DR (O 1 ^ 1 , id 3 ,^ 1 ). For inv + , define 

/ = DR (l ^ 1 , mul o jx mul o jx (id 1 ^ 1 , id 1 ^ 1 )) o id 2 ^ 1 ) , (10) 

inv+ = fo (addo jx (id 1 ^ 1 , -l 1 ^ 1 )) , (11) 

or, more colloquially, 

ft=l- /(/r) 2 dr, inv+t = /(t-l). (12) 

Jo 

Square root is defined analogously by 

ft=l+ [ inv + (2-/r)dr, sqrt + t = f(t - 1). (13) 

Jo 

Logarithm and exponentiation are analogous, using suitable integral equations. For the 
trigonometric functions, let sin = id 2 ,^ 1 o trig and cos = id 2 "* 1 o trig, where 

trig = DR (JX (O ^ 1 , l ^ 1 ), JX (id 3 "" 1 , (mul o jx (-l 3 ^ 1 , id 3 ^ 1 )))) , (14) 

which is to say, 

/sin A / 0\ f ( cost \ , 
The circle ratio is 7r( ) = 4- Arctanl, with Arctan defined by a suitable integral equation. □ 



Some authors say "the function 1/x is real primitive recursive" to mean that inv + is. It is 
not clear how such assertions without specification of domain can be justified. 

The reader may have felt uncomfortable with the unwieldy process of Definition [9] in picking 
the right solution h v out of H v . This can be simplified if we discuss only real primitive 
recursive functions, because of the following facts that result from the Uniqueness Theorem 
for initial value problems [H| and the Cauchy-Kowalewsky Theorem 10, Section 2.4]. 

Theorem 12. Let f , g, v, H v and h v be as in Definition® 

(a) If g is an analyti^ function with open domain, H v is the set of all restrictions of h v . 

(b) If f and g are analytic functions with open domain, so is DR(f,g). 

The fact (jgj) says that a solution of @ may diverge to infinity at some point but can never 
"branch" as in Figure HJ provided g is smooth enough. We therefore could have dispensed 
with Lemma fl~8l and simply let h v be the (graph) union of H v , so far as real primitive recursive 
functions are concerned, because they are analytic by (jbj). 



'This fact is often stated with a weaker assumption that g be Lipschitz continuous. 
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Integrate? 




Figure 2: Integrand with a singularity. 



3.3. Campagnolo's differential recursion 

The clauses of Definition guarantee that the integral equation © makes sense for all 

t G dom/i. The clause ([c]), however, could be slightly relaxed, since a small set of singularities 
in the integrand does not affect the integration. Define DRc by replacing ([c]) with 

([cj) {y, r, hr) G doing for any r G dom/i \ S, where S is a countable set of isolated points. 

This is due to Campagnolo jl, Definition 2.4.2], though he does not present a precise spec- 
ification of the "unique" solution as we noted in the Section 13.21 The choice between ([cj) 
and (jcj) is somewhat similar to the discussion regarding Principle [7J The issue is whether 
g(v,r,hr), where r G [0,t], is a "subexpression" of the right-hand side of the equation ©• 
Without going into the philosophical discussion to ask which is "natural," we point out some 
differences this choice incurs. 

Theorem [12] (jbj) fails if we replace DR by DRc, as the following example shows (Figure [2]). 
The function j:CR-*R defined by dom g = R\{1} and gt = inv + (sqrt + (sqrt + (t — l) 2 )) = 
1 /\/\t — 1| is real primitive recursive by LemmaHH But k = DR C (-2°^\go id 2 ^ 1 ) :C R -> 
R, where — 2 ^ 1 is the constant function with value —2, is the total function given by 



[-2- y/-t + 1 if t < 1, 

which is not different iable at 1. Note that dr(— 2 ^ 1 , g o idg - * ) is its restriction to (— oo, 1) 
and thus analytic. For a subtler example, recall the equation (1131) for sqrt + ; with DRc, the 
same equation produces the square root function defined on [0, oo), rather than on (0, oo). 

This breaks the assumption of Theorem [12] ([2$ and thus gives rise to incomparable functions 
in H v when, say, / = l ^ 1 and g = k o id 2 ^ 1 , with k from (TIBI) : that is, the equation 

ht = 1 + f k{hr)dT (17) 
Jo 

has two solutions that take different values at a point. 



Keeping the class analytic also conforms to Moore's intention [13|, Definition 9] to make 
the equation ([9]) equivalent to 

h0 = fv, 1 D 1 ht = g(v,t,ht), (18) 
which would not make sense for non-differentiable h. 



8 



3.4. A primitive recursive but not differentially algebraic function 

Claim [T] would not make sense if we adopted DRc in defining real primitive recursive functions, 
because there would then arise non-analytic functions, as we noted above. We now show that, 
even under our restrictive definition with the analyticity-preserving DR, the claim fails. 
Define f :C R 2 -> R by domT = (0, oo) 2 and 



f R 

r{R,x)= exp((x-l)-lnt-t)dt. (19) 

Jl/R 



*R 
l/R 

Define Euler's gamma function f:CR^Rby domT = (0, oo) and 



Fx = lim r(R,x). (20) 

It can be verified that this value converges and satisfies 

DTi= limD (M f(i?,i) (21) 

R^oo 

for each nGN and x G (0, oo). Holder showed that r is not differentially algebraic fy. 

We do not know if r is real primitive recursive, but t is easily shown real primitive 
recursive, using Lemma [TTJ However, contrary to Claim (TJ it is not differentially algebraic. 
For assume that it were. We would then have a nonzero polynomial P such that 

P(r(R,x),B^r(R,x),...,B ( -°' aiityP -^r(R,x)) =0 (22) 

for each (R,x) G (0, oo) 2 . Note that we used the characterization (JI|) of Theorem [3] in order 
to take P independent of R. We take the limit of (I22p as R — > oo, which by (1211) yields 

P(rx,D 1 rx,...,D arityP - 1 rx) = 0, (23) 

contradicting Holder. 

4. Other classes and related works 

This section discusses some other operators introduced by Moore and other authors. 

4.1. Minimization and Moore's real recursive functions 

For a function / :C R m+1 — > R, Moore defines MN / :C R m — >• R by 

Jf+ = inf{f>0 | /(«,«) = ()} ift+<-r, 

MN JV = < (24) 

It = sup {t < | f(v,t) — 0} otherwise. 

The class of real recursive functions^] is the smallest class containing all real primitive recursive 
functions and closed under jx, CM, DR and MN. 



5 This "recursiveness" of Moore's should not be confused with the same word also used in the context of 
Computable Analysis. As we see in AppendixEJ Moore's real recursive functions can even be discontinuous. 
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Moore states the definition of MN in a way that leaves ambiguous whether (124p has a value 
when, say, dom / = R m x [1, oo) and f(v,t) = 2 — t for all t > 1. Should it have the value 2, 
or be left undefined because "the zero-searching program gets stuck" ? 

It turns out that, whichever definition we choose, Moore's claim about iteration [l3|, Propo- 
sition 11] remains true, in the following modified form. Since the original proof again forgets 
partial functions, we present a new proof in Appendix O 

Lemma 13. If / :C R m — > R m is real recursive, there is a real recursive function g :C 
R m+1 — > R m that extends the function g' defined by dom g' = { (v , k) G R x (N \ {0}) | v G 
dom/ fc } and g'(v, k) = f k v for all {y , k) G doing', where f k = f°---°f- 

k 

We have to note, however, that the class of real recursive functions is probably not well- 
behaved, since, with MN producing non-smooth functions, the class no longer enjoys Theo- 
rem H21 We therefore doubt the significance of Claim [21 although it could be justified by 
using Lemma [13] to simulate Turing machines as Moore did. 



4.2. Linear differential recursion 

We have seen that many of the problems in Moore's original work were caused by failure to 
deal with partial functions properly. Some authors avoid this trouble by studying only oper- 
ators preserving totality, so that partial functions never come into discussion. Campagnolo 
and Moore take this path by considering linear differential recursion in place of DR. For 
classes defined by this operator, some relationships with digital computation are known 0,0]. 



4.3. Open problems 

Claim [1] has been the main rationale for calling variants of Moore's classes a model of analog 
computation. Now that we have lost it, an important challenge is the following. 

Open Problem. Find a subclass of our real primitive recursive functions, preferably with an 
equally simple definition, that has a close relationship to the differentially algebraic functions. 

Another direction would be to reformulate further the rest of Moore's work, as well as other 
authors' works that also suffer from the same kind of ambiguity. For example, it may be 



interesting to work out Mycka and Costa's class arising from the operator of taking limits [14] . 
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A. Old results 



We list some known theorems that we used in Section [2) 

The following Baire Category Theorem is used in the proof Theorem O 

Theorem 14. Let J be a subset of R m . The union of countably many closed subsets of J with 
empty interior has empty interior. 

Proof. Let Jo, Ji, . . .be closed subsets of J with empty interior, and U be any nonempty open 
subset of J. We will show that U \ UpeN ls nonempty For each P G N, we take xp £ R m 
and £p £ R as follows. Write B(x,e) for the open set of points in J whose distance from x is less 
than e. Let xq G U and eo £ (0 5 1) be such that B(xo,£o) ^ U. For each P £ N, let xp + \ G J7 
and ep+i G (0, 2 _p_1 ) be such that B(xp + i, £p+i) C B(xp,ep) \ Jp. This is possible because 
B(xp,£p) \ Jp is open and nonempty, since Jp is closed and has empty interior. As P tends to 
infinity, xp converges to a point in U \ UpeN Jp- ^ 

The proof of Theorem [3] also uses the following Identity Theorem (for real analytic functions 
of several variables), also known as the Principle of Analytic Continuation. It can be proved by 
straightforwardly generalizing the same assertion for unary functions 0, Section 1.2]. 

Theorem 15. An analytic function with open connected domain that vanishes on an open set 
vanishes everywhere. 

Let J C R be an open interval. It is well known that functions uq, . . . , Uk-i G C^fJ] are linearly 
dependent if and only if the determinant | (D l Uj)ij = o,... ! fe_i| , called their Wronskian, is zero. Using 
this fact, Ritt and Gourin [17|] showed (lurk) =>■ (lm|) of Theorem El 

Theorem 16. Let J CR be an open interval and let f G C W [J]. If we have 

P(/,D 1 /,D 2 /,...,D arit ^ p - 1 /) = (25) 

for some R- coefficient nonzero polynomial P , then we have (I25D for some Z- coefficient nonzero 
polynomial P. 

Proof. By the assumption, there is a finite set B C N antyP such that the functions 

f»° ■ (D/T 1 • • • (D arity p ~ 1 / ) ^ arity p_ 1 , for K . . . , ^it y p-i) G B, (26) 

are linearly dependent. The Wronskian of (|26|) thus vanishes, which is a Z-coefficient polynomial 
in /, D 1 /, . . . , D arity P +\ B \- 1 f\ This polynomial is nonzero, since otherwise (126ft would be linearly 
dependent for arbitrary /, which is absurd. □ 

One direction of Lemma [5] uses the following Transcendence Degree Theorem. 



Theorem 17. Let F be a subfield of a field E and D be a subset of E. If D C F(B) for some finite 
set B C E, then D C F{C) for some finite set C C D. 

Proof. For each d G D, the assumption gives 

d l =E(3 r d^ (27) 

j=0 
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for some I £ N \ {0} and (3j £ F(B). Suppose that for some d = do £ D, this equation contains 
some b £ B \ D, since otherwise we are done. Then we can rewrite (|27p as 

b^Zai-V (28) 

i=0 

for some k £ N \ {0} and a { £ F(B'), where B' = {B \ {&}) U {d }. 

For each d £ D and i G N, we can substitute ([271) and (i28j) repeatedly in d 1 to write 

= EE^r^' (29) 

i=0 j=0 

for some jij £ F(B'). The fc ■ Z + 1 elements 1, d, d 2 , . . . , d k l are hence linearly dependent over 
F(B'). We have thus found a set 5' with D C F(B') such that B' \ D has strictly less elements 
than B \D. Repeat. □ 

B. Maximal unique function 

This section shows that, from a set K of functions with a certain property, we can choose a function 
of which all functions in K is a restriction. This was used to justify Definition [9] in the presence of 
non-analytic functions where Theorem [T2] jaj) does not apply. 

We say that a set / C R is 0-convex if it is either the empty set or a possibly unbounded interval 
containing 0. Note that the union of 0-convex sets is 0-convex. 

We say that a set K of functions from R is consistent if for any t 6 R, the set {gt \ g £ K} 
has at most one element. In this case, the union of K means the unique function k such that 
dom/i = [J g£K domg and for each t £ dom/i, there is some g £ K with ht = gt. 

Lemma 18. Let H be a set of functions from R with 0-convex domain. Then the set K of functions 
unique in H is consistent. Moreover, if its union belongs to H , it belongs to K. 

Proof. For the first claim, suppose otherwise. Then there are functions £?o, &i G K and t £ 
domfco H domfci such that k§t ^ kit. This contradicts the fact that ko is unique in H. 

For the second claim, suppose that the union k of K is not unique in H. That is, there are a 
function g £ H and t £ dom/c n doing such that gt ^ kt. There is ko £ K for which t £ domfco- 
We have kot = kt ^ gt, contradicting the fact that ko is unique in H. □ 

This lemma can be applied to H = H v in the situation of Definition [9j because there the union 
of any consistent subset of H v belongs to H v . 

C. Iteration 

As we noted, the definition (|24|) of the operator MN is ambiguous, as it contains a subexpression 
f(v,t) that may be undefined for some (v,t). So when is MN fv defined? Possible answers include: 

(a) When t + and t~ are defined. 

(b) When at least either t + or t~ is defined; the condition t + < — 1~ will be used only when both 
are defined. 

And when is t + (resp. t~) defined? Possible answers include: 
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Figure 3: Simulating iteration f k v by real recursive functions. 

(i) When there is t > (resp. < 0) such that f(v,t) = and (v,t) G dom/ for all r G R. 

(ii) When there is t > (resp. < 0) such that f(v,t) = and (w,r) G dom/ for all r G [— 
(resp. [*, — *])- 

(iii) When there is i > (resp. < 0) such that f(v,t) = and (u,r) G dom / for all r G [0, i] 
(resp. [t,0]). 

(iv) When there is t > (resp. < 0) such that f(v,t) = 0. 

For (jl|), (|n]) and (|m|) . we may also consider adding the phrase "except for some countably many 
isolated r" (compare (jcf) in Section l3~3j) . 



Moore's informal explanation by a programming language 13|, Section 7] seems to suggest (jbj) 
and (jnj). However, without discussing which is the "right" definition of mn, we show that, whichever 
we choose, Lemma [T3l holds. The following proof is consistent with any of the above 2x7 possible 
definitions. 

Proof of LemmaMtA Denote MN fv by jxt. f (v,t). Let 

zero?x = /jy. (x 2 + y 2 ) ■ (1 — y), (30) 
integer? x = zero? (sin (tt ■ x)) , (31) 
rounds = x — fir. integer? (x — r), (32) 

so that ([32]) is the unique integer in (x — 1/2, x + 1/2]. We get inv of ((TJ) by 

imrx = fit. x ■ (x ■ t — 1). (33) 

The above four functions are total. Let 

digit (x, b, i) = round f -r ) — b ■ round f -— | (34) 

for b > 0, where b l = exp(i • In b). When b > 1 and i are integers, digit {x, b, i) is the digit in b l, s 
place when x is written in base-6 notation. Define 

elk* = digit (t, 2,-1), zigzagt = 0+ I (2 - 4 • clkr) dr, (35) 

Jo 

g(v,t)\ = (v\ + /V2 • (1 - clkr) • {f{h(v,r) - clkr ■ (h^) - v)) - h(v,r))\ ^ 
h(v,t)J \yj J \ 2 • clkr • (h(v,r) - g(v,r)) • inv(zigzagr) / ' 

as in Figured We have f k v = g(v, k - 1/2) for k G N \ {0}. □ 

Note that clkr • (h(v,r) — v) in ([36]) cannot be dropped, because of Principle [7] 
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